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ABSTRACT: Effects of plasma drift on the current collection by a long conducting 

cylinder in a magnetized plasma is studied by means of a 2 1/2 dimensional PIC code. It is 
found that for the drift velocity (V Q ) perpendicular to the magnetic field B Q , the electron 
current collected by a positive cylinder is considerably enhanced depending on the drift 
velocity. The distributions of plasma and the potential structure around the cylinder for 
several relative orientations between V Q and B are presented along with the comparisons of 
current with and without the magnetic field. Simulations with the magnetic field in the 
simulation plane show that the potential structures around the cylinder are 
two-dimensional double layers with dimension (L^) perpendicular to B much smaller than 
the dimension (Ly) parallel to B. In fact, L^ is found to be approximately determined by 
the current limiting radius given by the Parker— Murphy model. However, it is found that 
the collected currents in the simulations are generally higher than those given by this 
model. 

1. INTRODUCTION 

The knowledge of current collection by conducting bodies in space plasma is 
relevant to numerous applications such as the operation of plasma probes, charge 
neutralization on space vehicles, working of the solar cell arrays and the operation of an 
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electrodynamic tether. Most theories dealing with this topic are limited to simple 
geometries and idealized plasma models. For reviews of the theoretical studies, the reader 
is referred to Whipple and Laframboise and Sonmor in this volume. These reviews show 
that there is a general lack of theoretical studies on current collection in a magnetized 
plasma when there is a relative drift between the magnetized plasma and the current 
collector. The purpose of this paper is to contribute to this area by means of computer 
simulations using a PIC code. 

Our computer model is two dimensional; the axis of the cylinder is perpendicular to 
the plane of simulation. The magnetic field is oriented along the axis of the cylinder or in 
the simulation plane in different simulations, which bring out the effects of relative 
orientation between the magnetic field and the plasma drift on the sheath structure and the 
current collection properties. For the axial magnetic field, a simple-minded picture with 
radial electric fields indicates that the E * B drift will cause a magnetic insulation stopping 
any collection of electrons by the positive cylinder. Simulations show that this picture is 
not valid when there is a relative flow between the plasma and the cylinder; the potential 
structure is considerably modified so that the flowing electrons are focused onto the 
cylinder, making possible the collection of a relatively large electron current. 

When the magnetic field is in the simulation plane, the potential structure is 
extended along the magnetic field and its transverse dimension is quite limited and it is 
found to be given by the current limiting radius calculated by Parker and Murphy [1] in a 
non— flowing plasma. We find that when the flow is perpendicular to the magnetic field, 
the electrons intercepted by the extended field— aligned potential structure are partially 
collected by the cylinder and the current is found to be considerably enhanced over the 
current predicted by the Parker— Murphy model [1]. However, for the flow parallel to the 
magnetic field, the current is seen to be limited in a fashion described by the above model. 
Since in the low earth orbit, the orbital velocity vector is at large oblique angles with 
respect to the geomagnetic field, a current enhancement is expected. 
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2. SIMULATION TECHNIQUE 

The plasma flow past the cylinder is simulated as shown in Figure 1. The hatched 
area is the end view of the long conducting cylinder of radius r g . In the rest frame of the 
cylinder, plasma flows along the positive x direction with the velocity V Q . The flow is 
facilitated by imposing a dc convection electric field E q so that V Q = * B/B . The 

simulated plasma region is limited to r < R max (Figure 1). At the initial time t = 0, the 
simulation region is a vacuum and the plasma flow for t > 0 is maintained by injecting 
charged particles at the rim of the simulation box (r = R max ) over the angular region 7r/2 
< $ < 3^2 (Figure 1). The average injection velocity of the charged particles is V Q . At 
each time step a predetermined number (N-) of electron-ion pairs are injected to 
simulate a desired plasma flux. The plasma particles used in the simulations are like rods 
parallel to the axis of the cylinder [2]. The injected particles are chosen from Maxwellian 
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Figure 1. Geometry of the simulation. The conducting cylinder is shown by the 
hatched region. R gives the radius of the simulation system. Plasma 

IQSUC 

flows across the cylinder with a velocity V . The magnetic field is parallel to 
the axis of the cylinder. 
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distributions with electron temperature T = T and ion temperature T. = 0. The Y 

CO 1 

coordinates of the particles are chosen according to a uniform probability distribution. The 

2 2 X / 2 

X coordinates are first calculated by X = (R max - Y ) and then further randomized by 
replacing X by X + VAt where V is the particle velocity randomly chosen from a 
Maxwellian distribution and At is the time step. Our injection technique is quite similar 
to that described by Aldrich [6]. 

The magnitude of charge (q^), per unit length of such computer particles, is 
obtained by balancing the plasma flux into the simulation region at the injection boundary 
and the simulated flux due to the injection of the charge particles at each time step of 
duration At, giving 

llj = 2 R max eN o V o At / N inj (>) 

where N q is the ambient plasma density and e is the magnitude of the electron charge. 

The injection of equal numbers of electrons and ions insures that no net charge is injected 
into the system. 

The temporal and spatial evolutions of the plasma and fields are calculated by the 
self-consistent solutions of the equations of motions [2] of all the charged particles and the 
Poisson equation for the electric potential <fr. It is important to note that in our 
simulations, the electric field has two contributions as indicated by the following equation 

E = ^ + E x (2) 

where Eg is the convection field and E^ is determined by the space charges and the bias 
potential on the cylinder. Since Eg is uniform in space, V • Eg = 0 and the divergence of 
(2) gives 

V • E = V • Ej = p/e Q (3) 

where p is the electric charge density. Under the electrostatic approximation, E^ = — V0 
and (3) gives the Poisson equation 

V 2 * = -P/' 0 ( 4 ) 
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The boundary conditions on the electric potential <j> are <p(r = r g , 0) = (f> Q and <p ( r = 
R max’ ^ = °’ where is the bias P° tential of tbe c y Under - The Particles striking the 
cylinder and those leaving the system are assumed to be lost. However, the simulation 
system is maintained quasi— neutral at the ’global’ scale. For this purpose, we compare the 
total numbers of electrons and ions in the entire system at each time step. The deficit 
charged particles, which are taken from a Maxwellian plasma reservoir are randomly 
distributed over the entire simulation system according to a uniform probability 


distribution. 

The collected current (I) is calculated by counting the electrons and ions striking 
the cylinder during each time step, 


1 = s a «y At (*> 

where is their number, and is given by (1). We note that although q Q depends on 

the numerical factors R , N- . and At; the current I is found to be independent of them, 

max’ mj 

if R and N- • are sufficiently large and At is sufficiently small. This was verified by 
max inj 

carrying out simulations by varying these parameters. 


3. NORMALIZATIONS AND DEFINITIONS 

We discussed earlier that the charge on a computer particle is given by (1). If 
qje = tj, the analogy between the real and computer particles requires that the masses m, 
effective temperatures T and density N satisfy the relations 

m c a = "“ra' T ca = ^ra and N ca =N tJ t > 3 * * (6) 

where the subscripts r and c refer to the real and computer particles, respectively. It is 

worth mentioning that the electron and ion Debye lengths and plasma frequencies are 

invariant under the scaling law described by (6) [4]. 

The results presented in this paper are based on simulations with the following 

1 1 —3 

ionospheric plasma parameters: ambient plasma density N Q = 10 m , electron 
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temperature T g = 0.2 eV, plasma debye length ~ 1 cm, electron plasma frequency Wp Q 

7 

1.8 * 10 rad/s, and the magnetic field B q = 0.3 G. With the above ambient plasma 


parameters, the electron thermal current J = N eV, /V ~ 1.2 mA/m, where V, = 

r O X6 16 

(kT 0 /n. e ) 1/2 ~ 192 km/s. 

In order to simplify the equations and to generalize the applicability of their 
solutions to different situations with varying plasma and current— collector parameters, we 


use the following normalizations: potential (j) = <P n = kgT e /e; time t = tu/p Q ; 

velocity V = V/V^ e and distance r = r/A^. 

In view of the above normalizations the Poisson equation (4) can be written as 


d 2 ~<t> 1 d 2 (p 1 d 2 '<j) 


T? 


+ -■ 


+ 


r di r d<t>‘ 


'do 


<f>. 


(q n ic - <l n prJ 


ec 7 


( 7 ) 


n 


where qn^ and qn gc are the charge per unit volume associated with the computer ions and 

electrons, respectively. It is assumed that both types of particles have the same magnitude 

of charge, i.e., q = q. = q, as given by (1). These charge densities (qn. and qn ) are 

determined by calculating the number of computer particles at each grid point by the area 

sharing method [2] and dividing it by the effective volume of a cell. This volume is given 

by r-A0ArAz, where r- is the radial distance of a grid point, Ar and A0 are the radial and 
J J 

angular grid spacings, respectively and Az is the length along the axial direction. With 
these definitions and equation (1), the normalized Poisson equation takes the form 


^f + l.^i+'^L = _2E 




+ f “ + ? W* ■ ■ 2R m ax V o A ‘( Al 'ic- A “ec)/ r j ArA ®inj 


( 8 ) 


where An- and An are the number of computer ions and electrons shared on a grid point 
1C 6C 

(jAr, iA0). The above equation is solved by employing FFT in 6 and triadiagonal method 


in r. 


The numerical results presented below are based on the following numerical 
parameters: R- max = 140A ( j » 1.4 m, r g = 10A^ o s 0.1 m, At = 0.2, Ar = 1, A0 = 10‘ 
and the normalized flow velocity V = V /V. is varied. The simulations are carried out 
with H + plasma for which nij/m e = 1836. We note that in our simulation electron 
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cyclotron period r cg s 18 while the ion cyclotron period - (m./m e ) t cg . Thus, for 
the time scales in the simulations, electrons are magnetized, while ions behave as 
unmagnetized charged particles. 

4. NUMERICAL RESULTS 

In the following discussion we first present results for B = 0, which are later used 
for the purpose of comparisons with the results for non — zero magnetic fields with different 
orientations. 

4.1 Simulations with B = 0 

We recall that the simulation begins with no initial plasma in the system. The 
simulation plasma builds up in the system in response to the injection of particles as 
described above in Section 2. In the simulation described in this subsection <p Q = 100 and 
V Q = 0.3. Figures 2 and 3 show the evolution of the plasma; in Figure 2 the contours of 
constant density of ions are shown at some selected times. The minimum density contour 
is n = 0.1 and the density interval between the contours is An = 0.3. The electron density 
shows nearly the same evolution as the ions. After about t = 550, a quasi— steady state is 
reached in the plasma distribution. 

The distribution of the computer particles in the r — 0 plane are shown in Figure 3, 
each dot in the panels of this figure represents a particle. The left-hand panels show 
electrons while the right-hand ones show ions. Note the formation of a distinct wake 
behind the cylinder (also see Figure 2). Another noteworthy feature of Figures 2 and 3 is 
that a bow structure forms in the ram direction; in this structure the density is generally 
enhanced. We also see from Figure 3 that ions are not able to reach the cylinder because 
the kinetic energy of the ions (1/2 m- = 82.6 kT Q ) associated with the plasma drift is 
smaller than the potential energy e<p Q = 100 kT Q , where <f> Q is the bias voltage on the 
cylinder. Outside the wake region, the plasma density n = 1, indicating a uniform plasma 
flow in the ram direction away from the bow structure. 
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Figure 3. 


Evolution of plasma is shown by showing the distribution of (a) electrons and 
(b) ions in the r-0 plane. Each dot in this figure represents a computer 

particle. <p Q = 100, B = 0, V Q = 0.3. 
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The evolution of the potential distribution around the cylinder is shown in Figure 4, 
which gives the equi potential surfaces at some selected times. The contour levels are A0 = 
5 apart. This figure shows that after about t = 500, the sheath structure reaches a 
quasi-steady state. In the wake region the potential is generally negative. 

In response to the evolution of the plasma and potential around the cylinder, the 
collected current evolves as shown in Figure 5. The current reaches a quasi— steady state 
after about t = 500, in agreement with the evolutions of the density and the plasma 
potential. After this time, the plasma and the potential are still undergoing some changes, 
especially in the wake region, but they seem to have negligible effect on the current 
collection. The time— average current for B = 0 in the quasi— steady state (t > 500) is 
about I ~ 18 mA. 



X 

Figure 4. Evolution of the equipotential surfaces around the cylinder. Note that after 


about t = 600, the equipotential surfaces attain a quasi— steady state. The 
equipotential contours are = 5 apart. = 100, B = 0, V Q = 0.3 
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TIME (tOpe) 

Figure 5. Temporal evolution of the current collected by the cylinder. The thick line 
curve shows the time average current when fast oscillations are averaged out. 
Note that the current attains a quasi— steady state when about t > 600. 
'<t> 0 = 100, B = 0, V Q = 0.3 

The simulations with B = 0 were carried out for several bias potentials. Figure 6 
shows the V— I characteristics of the cylinder. It is found that I a which is in 

agreement with the orbit— limited current collected by a cylinder (e.g. see Chen [5]). 
However, the proportionality constant is found to be given by 

I ~ 1.8 (e^/kTg) 1 / 2 mA/m , 

which is found to be by a factor of two larger than that for V Q = 0. It is expected that in 
the limit V 0 = 0, the simulations must yield the current as predicted by the orbit-limited 
current. However, the simulation runs with very small drift velocities take too long to 
complete and so far we have not carried out such simulations. 
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Figure 6. I-V characteristic of the cylinder. Note that the horizontal axis is <p Q ' 
The current variation shows the linear relation I a = 0.3, B = 0. 


4.2 Simulation with B = B z = 0.3 G and = 100 

We do not show the temporal evolution of the plasma and potential here, instead we 
just present here the quasi— steady state distribution of the plasma and potential around 
the cylinder. The top two panels of Figure 7 show the contours of constant ion and 
electron densities. The corresponding distributions of the particles in the r — 9 plane are 
shown by the two middle panels. The bottom single panel shows the distribution of 
potential; equipotential surfaces at intervals of A0 = 5 are shown. It is worth pointing out 
that the plasma and potential distributions for B z = 0.3 Gauss is quite different from those 
for B = 0. In the former case (B = 0.3 Gauss), the equipotentials show a multicell 
convection pattern [6]. The fan— shaped structure extending below the cylinder is the 
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consequence of the stagnation of the plasma flow below the cylinder due to the opposition 
to the flow by the E * B drift in the initial radial electric field [6]. The fan— shaped 
equipotentials cause electrons to circulate around the cylinder due to the E»B drift. The 
electron flow coming from the left is caught in this convection cell and focused on to the 
cylinder as shown by the crowded equipotentials immediately on the top of the cylinder. 
This circulation of the flowing electrons facilitates their collection by the cylinder. 

The temporal evolution of the current collected by the cylinder with the axial field 
B is shown in Figure 8. The current is seen to reach a quasi-steady state at about 
t ~ 700, after which its average value I ~ 14 mA/m, which is only slightly lower than 18 
mA/m for B = 0. The simulation with the axial magnetic field shows that the magnetic 
insulation due to E * B drift in the initial radial electric field is destroyed due to the 
considerable modification of the potential distribution caused by the plasma flow [6]. 



TIME (CDpe) 

Figure 8. Temporal evolution of the current for B z = 0.3 G, <f> Q = 100 and V Q = 0.3. 
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Simulation with B = By — 0.3 Gauss, and 50 

Note that in this case the magnetic field is in the plane of the simulation. This 
allows us to study the B field-aligned potential structures. The quasi-steady state 
distributions of the plasma and potentials are shown in Figure 9. The top two panels show 
the contours of the ion and electron densities. The corresponding distributions of ions and 
electrons in the r-0 plane are shown in the middle panels. The wake structure is clearly 
seen from these panels. In the ram direction the plasma is generally uniform with the 
normalized density n = 1. The bottom panels show the equipotential surfaces from 
simulation with 4> 0 = 50 (left) and \ = 25 (right). These bottom panels show that the 
potential distributions are extended along the magnetic field. When = 50, the potential 
structure is seen to extend all the way to the boundary of the simulation plasma. In order 
to examine the effect of the boundary the simulation was repeated by lowering </> 0 to 25 
and increasing the size of the system from R max = 140 to 185. The result is shown in the 
bottom right-hand panel. It is seen that potential structure is now nearly fully 
accommodated in the simulation region. 

It is interesting to examine the size (L^) of the potential structure transverse to the 

magnetic field. Figure 10 shows the radial distribution of the potentials for 4> Q = 50 and 25 

in the ram direction (0 = 180*). It is seen that the potential structure becomes narrower 
with increasing <f> . The radial distances at which </> = 0 for the above bias voltages are 

given by 

L ± ~ 25 A d , 4> Q = 25 ( 9a ) 

L x ~ 31 A d , 4> q = 50 (9b) 

Parker and Murphy [1] have considered the collection of electrons by a positive 
sphere. Using conservation of energy and angular momentum, they have shown that in the 
case of non— flowing plasma, the electrons which are possibly collected by the sphere, are 
confined in a cylinder of radius r Q as shown in Figure 11, where r Q is given by 





Quasi— steady state distributions of (a) electron density, (b) ion density, (c) 
electrons, (d) ions, (e) potential for = 50, V Q = 0.3, B = B^. = 0.3 Gauss, 
and (f) potential distribution for = 25, V Q = 0.3, B = B^. = 0.3 Gauss in 
a simulation with larger system size. 








where p Q( p is the electron Larmor radius with the electron energy e<p Q . It is interesting to 
note that the value of L estimated above for ~<t> Q = 25 and 50 are quite accurately given by 
(10), with a as the radius of the cylinder. 



o 40 80 120 
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Figure 10 Radial distributions of potential for <f> Q — 50 and 25 in the ram direction 
(6 = 180°). B = B = 0.3 Gauss. V = 0.3. 





Undisturbed Plasma 

Figure 11 Geometry of the Parker— Murphy model for the current limitation. Electrons 
contained inside the cylindrical volume of radius V Q are possibly collected by 
the electrode at a positive potential. 


According to the Parker— Murphy model the current collected by the spherical 
electrode in a non— drifting plasma is simply the electron flux intercepted by the cylinder of 
radius r Q (Figure 11): 


hit = J r 


( 11 ) 


where J r is the electron thermal current given by J r = N q eV te /V In the present 
situation, is associated with the thermal motion along the y direction parallel to B . 

y 

We find that for the flow in a direction transverse to the extended potential 
structure, the current is considerably enhanced. If the current was collected primarily 
through the two ends of the potential structures (Figures 9e and 9f), the total electron 
current collected by the cylinder is given by 


*PM _ ^ x 2r Q J f A/m (12) 

which is only about 1.2 mA/m for = 25. Our simulation shows a considerably larger 
collection of electron current. Figure 12 shows the evolution of the current collected by the 
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cylinder when 4> Q = 25; in the quasi-steady state the current is about 7.5 mA/m, which is 
found to be close to the current collected without the ambient magnetic field with the same 
drift velocity V Q = 0.3 (Figure 6). The excess current (~ 6.3 mA/m) is interpreted in 
terms of the interception of the electron flow by the extended potential structure along the 

magnetic field. 



4.4 Simulation with B = a x B Q 

In order to examine the effects of relative orientation of the drift velocity with 
respect to the ambient magnetic field B^ in the simulation plane, we carried out another 
simulation with V Q || B Q . The potential structure for this case in the quasi-steady state is 
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shown in Figure 13. It is seen that now the potential structure is extended along x, the 
direction of the B field. Its transverse dimension is again found to be limited according 
to (10), which gives the current limiting Parker— Murphy radius as a function of the bias 
voltage <f> Q . The temporal evolution of the collected current for B = a x B q is shown in 
Figure 14a. For the purpose of comparison, the current with B y is plotted in Figure 14b 
for the same value of ^ = 50. Note the different vertical scales in Figures 14a and 14b. It 
is seen that for the flow along B, the current is significantly reduced compared to the case 
with flow transverse to B. As noted earlier, in the later case the interception of the flow by 
the elongated potential structure enhances the current. 

It is instructive to quantitatively compare the current from the Parker-Murphy 
model with that from the simulation with B x> We already saw that the former current is 



Figure 13 Quasi— steady state potential distribution for B = B x = 0.3 Gauss, V Q = 0.3, 

h = 50. 
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given by (12). However, a few observations must be made while using this expression in 
the present case. The flow along the magnetic field makes the potential structure 
asymmetric with respect to x = 0 because of the formation of the wake behind the cylinder, 
making the first factor of 2 in (12) inaccurate. The use of this factor will give an 
overestimate. Furthermore, J r in (12) must be replaced by a modified current density due 
to the plasma flow; for the flow velocity V Q = 0.3 V tg , this modified current density 
J g = 1.44 J f . With these considerations, (12) gives I pM < 2.2 mA/m. Figure 14 shows 
that the time average current is about 3.5 mA/m, which is at least 60% larger than the 
current predicted from the Parker— Murphy model. The enhancement in the current 
suggests the transport of electrons across the magnetic field line. The exact mechanism for 
the cross— field electron transport has not been identified from the simulations. However, 
cross— field diffusion due to the fluctuations in the field need to be examined [7]. 

5. CONCLUSIONS AND DISCUSSION 

The main conclusions of this paper are as follows: 

(i) When the relative plasma flow is transverse to the magnetic field, the current 
collected by a positive electrode can be considerably enhanced depending on the 
relative drift velocity. 

(ii) For the flow along the magnetic field, the current is limited as predicted by the 
Parker— Murphy model [1]. 

(iii) Simulations with B in the simulation plane show that the magnetic 
field— aligned potential structure is like a double layer with dimensions transverse to 
the magnetic field determined by the limiting radius given by the Parker-Murphy 
model [1], but it is extended along the field line. 

(iv) Simulation with the axial magnetic field shows that the potential structure 
represents a multi— cell convection pattern. The fan— shaped cell is seen to focus the 
electron flow on to the cylinder thus destroying the magnetic insulation effect. 
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Figure 14 Temporal evolution of the 
h = 50. (b) B = B = 0.3 


irrent for (a) B 


Jauss, V = 0.3, 





(v) In the low earth orbit, the orbital motion is nearly perpendicular to the 
magnetic field, and the relative flow velocity V Q 2 8 km/s, which can considerably 
enhance the current collection. 

* In a real situation with an arbitrary orientation between ^ and V Q , the potential 
structure will be the combinations of the structures shown in this paper. 

Recently Myers et al. [8] have demonstrated that the measured currents in a rocket 
experiment agree well with the predictions from the Parker— Murphy model [1]. Raitt et. al 
[this volume] have carried out a similar comparison. Since in the rocket experiments the 
relative drift velocities are only a few hundred meters per second, the current enhancement 
due to the relative drift is not expected to be significant. 

In our present simulations plasma flows while the electrode is standing. In space 
the electrode cuts across the magnetic field lines. This raises some question about the 
dynamical effects. We note that in the simulations starting with an initial vacuum state, 
the quasi-equilibrium is reached quite quickly in a time of about 500 2 30 jj&. In real 
situations of space the quasi— equilibrium condition are expected to reach in a considerably 
shorter time. On the other hand, the contact time of a current collector with a magnetic 
flux tube depends on its sheath size. If we use the sheath size as given by (10) for large 
electrode voltages, the contact time can be estimated by 

r c = 2 ' r V / V o 

Using typical parameters (<p Q = 100 V, a = 1 m, B = 0.3 G and V Q = 8 km/s) it is 
found that t q > 250 fjs. Comparing this time with the sheath establishment time of the 
order of a few tens of microseconds, it is inferred that the quasi-equilibrium conditions for 
the potential structure and the current collection as found from the simulations are likely 
to be maintained for a current collecting electrode in the low earth orbit. 
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